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INTRODUCTION 

The purpose of this study was to select weather variables that could be 
used to estimate the soybean yields for the country of Brazil. Soybean produc- 
tion in Brazil had its inception in the state of Rio Grande do Sul where 
soybeans were first planted in rotation with wheat. Presently, soybeans are 
planted in seven southern states: Rio Grande do Sul, Santa Catarina, Parana, 

Sao Paulo, Mato Grosso, Goias and Minas Gerais. This means that a variety of 
growing conditions are covered in this extensive area. Figure 1 is a map of 
the soybean-growing areas of Brazil. 

The northern part of the soybean area, including the states of Mato Grosso, 
Goias, and Minas Gerais, are characterized by a "wet-and-dry" climate with a dry 
season of two or more months. The rest of the area is subtropical humid with 
relatively abundant rainfall which is well-distributed throughout the year but 
with slightly more rainfall in the warm months. Summers are hot and winters 
relatively mild. Temperatures in the summer are above 40°C in the plains of 
Rio Grande do Sul when warm air masses penetrate the plains. Southern Sao 
Paulo is the northern limit for frost occurrence. 

The soybean growing season begins with planting in October and November and 
runs through April and May when harvest starts. 

METHOD 

Multiple regression analysis of yield with selected agroclimatic indices 
was used tc derive a suitable model. The first approach taken was to develop a 
model for each state. However, lack of data for several states posed a problem 
and no suitable models were derived. Therefore, it was decided to create a 
covariance model whereby all the available data are combined to obtain a model 
for the country. 
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Figure 1. Soybean-growing areas of Brazil 
(Pitcher, 1971) 









Since moisture stress was considered a prime determinant of yield, the 
index P-PET (precipitation minus potential evapotranspiration) was used in the 
regressions. 

The regression equation is: 

A 

Y * oc+ BjD + B2R0FNi + B3<P-PET) i + E 

where 

A 

Y = Estimated yield, 

CC = Constant, 

Bj * Coefficients of the variables j * 1 - 3, 

D = Dummy variable to adjust each state's yield to a base yield set by 
Rio Grande's do Sul's yield, 

RDFNi - Deviation from normal of total precipitation for month i, 

P-PETi * Precipitation minus PET for month i, and 

E * Unexplained error 

In developing the model, various procedures of the Statistical Analysis 
System (SAS Institute Inc., 1979) were used. The procedures used and the 
operations performed with each are summarized in the Appendix. The .elected 
model had the highest R^ and included variables significant at (or close to) 
the 10 per cent level and agronomical ly meaningful. 

DATA 

The Brazil crop data for 1961 thru 1977 were obtained from the Foreign 
Agricultural Service (Sam Ruff, personal communication, 1982). The data were 
recorded with year of yield as year of harvest, so the weather influencing 
the crop occurred during year-1. 

Meteorological data from 1975 through 1980 were used to model because 
there is no apparent trend in the yield data during this period. Table 1 lists 
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• BrsgHle/Cruzelro 
• Golanfe • Montes Cleros 



Figure 2. Location of Meteorological Stations Used to Derive 
Data Sets for the Brazil Soybean Model. 
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STATE 

METEOROLOGICAL STATION 

WMO NUMBER 

Goias 

Brasilia/Cruzeiro 

83377 


Goiania 

83423 

Mato Grosso 

Ponta Pora 

83702 - 


Pres Prudente 

83716 


Foz do Iguacu 

83826 

Minas Gerais 

Montes Claros 

83437 


Araxa 

83579 

Parana 

Londrina 

83766 


Santa Branca 

83810 


Foz do Iguacu 

83826 


Porto Uniao 

83864 

Rio Grande do Sul 

Ijui 

83805 


Veranopolis 

83806 


Farroupilha 

83807 


Sao Borja 

83901 


Pas so Fundo 

83914 


Julio de Castilhas 

83935 


Santa Maria 

83936 

Santa Catarina 

Porto Uniao 

83864 


Sao Joaquim 

83920 


Vacaria 

83918 

Sao Paulo 

Pindorama 

83664 


Ribeiro Preto 

83668 


Mococa 

83680 


Bauru 

83722 


Jau 

83723 


Limeira 

83728 


Campinas 

83729 


Monte Alegre do Sul 

83730 


Tiete 

83777 


Sao Parlo 

83800 


Jundiqi 

83802 
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Table 1 


Meteorological Stations Used to Derive Data Sets 
for the Brazil Soybean Model. 


the stations used to derive the meteorological data set for each state. 

Figure 2 shows the location of each station. The seven states of Rio Grande 
do Sul, Santa Catarina, Parana, Sao Paulo, Mato Grosso, Goias, and Minas Gerais 
were included in the covariance model. 

PROCEDURES 

. The original variables used in the regression equation included "dummy 
variables" for six of the seven states (Rio Grande do Sul was the exceptir, 

The "dummy variables" adjust the contributions to yield of each of the states 
to a base yield which, in this case, is Rio Grande do Sul's yield. Weather 
variables selected reflected available moisture for the months December through 
March. The variables which were significant at the 10 per cent level were the 
"dummy variables" for Sao Paulo, Parana, and Santa Catarina and the P-PET for 
January and February (see the Appendix for a definition of P-PET). The coef- 
ficient for the Santa Catarina "dummy variable" was negative, indicating its 
yield was below that of the norm set by Rio Grande do Sul; the coefficients 
for the other dummy variables were positive. The coefficients for January and 
February P-PET were positive, indicating a need for moisture during those 
critical months. The statistics of the selected model are summarized in Table 2 

TEST RESULTS 

A jackknife test was run on the final model. In this test, a year was 
eliminated from the crop data and the model was used to predict that year's 
yield. This process was done for each successive year beginning with 1975. 

The test had to be run separately on each state. The results are printed on 
Tables 3 through 9 and plotted on Figures 3 through 9. 



APPENDIX 


Definition of Variables 

P-PET, precipitation minus potential evapotranspiration, is an index used 
to measure the amount of moisture available for plant growth. Potential 
evapotranspiration is determined by the procedure developed by Thornthwaite 
(1948) which uses only temperature: 



where I * Heat index, which is the sum of the 12 monthly indices i. 



T = Monthly temperature in °C, and 

a * An empirical exponent, 6.75 x 10"^I^ - 7.71 x 10"^i2 + 1.79 x 10~2 I 
+ 0.49. 

The duration of daylight is used to adjust potential evapotranspiration as a 
portion of 12 hours. 


Statistical Analysis System Procedures Used 


PROC CORR 

PROC PLOT 
PROC STEPWISE 


PROC STEPWISE FORWARD 


Computes correlation coefficients between 
variables, including Pearson product-moment 
and weighted proc.ct-moment correlation. 

Graphs one variable against another, producing 
a printer plot. 

Provides five methods for stepwise regression. 
Stepwise Is useful when selecting variables to 
be Included in a regression model from a 
collection of Independent variables. 

Segins by finding the one-variable model that 
produces the highest R^. For each of the other 
independent variables, FORWARD calculates F- 
statistlcs reflecting the contribution to the 
model if the variable were to be included. 


PROC STEPWISE BACKWARD Begins by calcu’ating statistics for a model 

including all the independent variables. The 
variables are deleted from the model one by 
one until all the remaining variables produce 
F-statistics significant at the .10 level. 


2 . . 


7 



PROC STEPWISE STEPWISE The stepwise method Is a modification of the 

forward selection technique, differing in 
that variables already in the model do not 
necessarily stay there. After a variable is 
added (as in the forward selection method) the 
stepwise method looks at all tne variables 
already included in the moael and deletes any 
variable that does not produce an F-statistic. 
significant at the .10 level. Only after this 
check is made and the necessary deletions 
accomplished can another variable be added to 
the model. 

PRUC STEPWISE MAXR (Maximum R^ imp nvement) Unlike the three 

techniques above, this method does not settle 
on a single method. Instead it looks for the 
“best" two-variable model, the "best" three 
variable model, and so forth. 

PROC PETM Uses latitude and mean monthly temperature to 

calculate Thornthwaite's potential evapo- 
transplration for each month. 

PROC ZINDEX Uses monthly PET’s, precipitation, SS 

(begi ning moisture in surface layer), 

AWCS (available water capacity in surface 
layer), SU (beginning moisture in the 
underlying layer), and AWCU (available water 
capacity in the underlying layer) to calculate 
Palmer's sojQ moisture budget, drought index 
Z, ET, and ET. 
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Table 2. Statistics of Brazil Soybean Model 
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Figure 4. State of Mato Grosso Soybean Model. 
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Figure 5. State of Minas Gerais Soybean Model. 
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Table 6. Results of Jackknife Test for State of Parana Soybean Model. 
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Table 7. Results of Jackknife Test for State of Rio C’.ande do Sul Soybean Model. 
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Figure 7. State of Rio Grande do Sul Soybean Model 
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Santa Catarina Soybean Model. 
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Figure 8. State of Santa Catarina Soybean Model. 
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Table 9. Results of Jackknife Test for State of Sao Paulo Soybean Model 
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Figure 9. State of Sao Paulo Soybean Model. 



